Ab initio investigation on oxygen defect clusters in UO2+ 
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By first-principles LSDA+U calculations, we revealed that the current physical picture of defective 
uranium dioxide suggested solely by neutron diffraction analysis is unsatisfactory. An understanding 
based on quantum theory has been established as a thermodynamical competition among point 
defects and cuboctahedral cluster, which naturally interprets the puzzled origin of the asymmetric O 
and O interstitials. It also gives a clear and consistent agreement with most available experimental 
data. Unfortunately, the observed high occupation of O site cannot be accounted for in this picture 
and is still a challenge for theoretical simulations. 

PACS numbers: 61.72.J-, 71.15.Nc, 71.27. +a 
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Fluorite structure commonly appears in rare earth and 
actinide oxides. One of them UO2 is the nowadays widely 
adopted nuclear fuel and most of whose applications are 
highly related to atomic defect behaviors. The current 
physical picture of defective UO2 is based purely on neu- 
tron diffraction measurements. It suggested that no oxy- 
gen can occupy the octahedral site (Oi) of the cation 
FCC lattice and all oxygen interstitials should displace 
from this site along the (110) and (111) directions about 
1 A(sites and ), respectively, to form the so-called 
Willis type clusters by associating with the nearby oxy- 
gen vacancies (0„).[l[ H, H, 0] Note the clustering inter- 
pretation in this picture is suggestive and might be am- 
biguous since explicit atomic positions are unavailable in 
those experiments. But this suggestion invalidated the 
simple point defect model and diffusion mechanism that 
widely employed to describe fuel behaviors. 0, H Also 
itself is inconsistent in that it requires a quite different 
cuboctahedral (COT) cluster to account for the closely 
related U 4 9 /U 3 07 phase. i, ©, El El El This sit- 
uation provoked confusion because it is difficult to in- 
vestigate the Willis clusters thoroughly and meanwhile 
people are unsure how large the clustering effects should 

be. 0, El, El El M 

On the other hand, as a powerful quantum me- 
chanics technology, density functional theory provides 
an accurate ab initio method to understand material 
behaviors. [Tt} For example the bulk stoichiometric UO2 
has been well described by local de nsity approximation 
with Hubbard correction (LSDA+U (Hf) functional, [ll] 
Application of this method to non-stoichiometric UO2+X 
also showed amazing agreement with available experi- 
ments, namely, in (i) negative volume change induced by 
oxygen interstitials, (ii) predominance of oxygen defects 
when x > 0, (hi) exclusive COT clusters at high x region 
(U4O9/U3O7) . 0, i, 1 1 1, S 0, S S EI HI In this let- 
ter we will investigate the full temperature-composition 
region to study the oxygen defect clustering behaviors in 
UO2+Z with this method. 



All defects were modeled in a 2 x 2 x 2 fluorite cu- 
bic cell with otherwise 96 atoms. Larger supercells also 
tried to check the size effect. Structure optimization has 
been performed with PI symmetry and until all forces 
and stress less than 0.01 eV/A. The employed compu- 
tational technique and parameters, such as the kinetic 
energy cutoff and sampling k-points, are the same as in 
Ref.0. 
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FIG. 1: Relaxation process of the Willis 2:2:2 cluster (a) to a 
V-30 split interstitial (b): each cube indicates one oxygen 
cage and the solid arrows point to atomic relaxation direc- 
tions; in (a) the topside interstitial will annihilate the nearest 
vacancy eventually and the two O interstitials move down- 
wards to the nearby O sites. 

Willis-type cluster has a geometry of that O in- 
terstitials pushing the nearest lattice oxygen displaced 
along the (111) direction and thus created a connected 
0':0„:0" defect cluster. 0, A typical 2:2:2 cluster is 
shown in figure QJa). Other already suggested clusters 
include 1:2:2, 4:3:2, 2:3:2 and so on.[13| We have checked 
the structural stability of several simple clustering mod- 
els in Ref. 0, but failed to give a definite conclusion for 
Willis clusters due to the limitation of the method em- 
ployed there. Here we will investigate the stability of 
1:2:2, 2:2:2 and 3:3:2 (with 2 inequivalent configurations, 
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TABLE I: First principles results for oxygen defect clusters 
in UO2+1: AV the defect induced volume change per fluorite 
cubic cell; E/, E e / and E p / are the overall defect formation 
energy, formation energy per excess oxygen and per oxygen 
Frenkel pair, respectively. 
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FIG. 2: Defect concentrations of point oxygen interstitial, 
uranium vacancy, COT-o, V-30 and V-40 clusters at 
1500 K, respectively. All other components are negligible. 



respectively), 4:3:2, and 2:3:2 (with 3 inequivalent config- 
urations) clusters directly. The behaviors of large Willis 
clusters can be inferred from these small ones. Calcula- 
tions showed that for each kind of cluster, all inequivalent 
configurations always relaxed to the same final structure. 
In details, all 1:2:2 clusters relaxed to a point oxygen 
interstitial occupying the octahedral site (Oi), all 2:2:2 
and 3:3:2 clusters relaxed to a V-30 cluster, the 4:3:2 
relaxed to a V-40 cluster, and all 2:3:2 decayed to point 
Oi. There is no energy trap presented along the relax- 
ation paths, showing that Willis-type cluster is not even 
met ast able. 

Figure[lja) shows the relaxation process of a 2:2:2 clus- 
ter, in which the arrows indicate the atomic movement 
directions. Figure [TJb) gives the final V-30 cluster. 
The failure of Willis-type clusters is owing to that O in- 
terstitial in fact cannot push the nearest lattice oxygen 
out along the (111) direction. Instead itself will relax to 
an O" site to form a split-interstitial (V r -3O"/^ / -4O")0| 
or to annihilate nearby vacancy. This is different from 
that in COT cluster where the integral support stabilizes 
the O interstitials. [2(| In the case of 2:2:2, the topside 
O will annihilate the nearest 0„, with a little distor- 
tion along the (001) direction; the two O s will move 
downwards to the nearest O sites and form a V-30 
cluster with the low O interstitial. This rule applies to 
other Willis clusters too. For example, In 1:2:2 and 2:3:2 
cases, O interstitials were not capable of supporting the 
structure and all interstitials annihilated nearby vacan- 
cies with the last one relaxed to an Oi site; the four O s 
in 4:3:2 also couldnot prevent the O s from annihilating 
and themselves then relaxed to O positions surrounding 
the middle vacancy and formed a V-40 cluster. 

This makes the concept of Willis cluster generally 
failed. For clusters bigger than 4:3:2, there is no long- 
ranged interaction that can stabilize the structure. We 
have tried several primary calculations on big clusters 
but failed to find any evidence that they were favored. 

Since all Willis- type clusters will relax to Oi, V-30 , 
V-40 and their combinations, it is reasonable to assume 
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that V-30 and V-40 clusters might take the role of 
Willis clusters to account for the measured occupation 
number of O and O sites. 0, [H, H| To verify this as- 
sumption, we calculated these two clusters in a 2 x 2 x 2 
supercell. The defect induced volume changes and en- 
ergetic information arc listed in table |TJ Like Oi and 
COT clusters, V-30 and V-40 clusters all result in a 
negative volume change, but is much smaller than that 
of COT-o (COT with one O" occupies the center). @, [1| 
As suggested by counting the nearest vacancy-interstitial 
pairs, [7( V-40 has lower overall and per excess oxygen 
formation energies. But due to the number of Frenkel 
pair available when compensating with point 0^, the for- 
mation energy per Frenkel pair of V-40 is a little higher 
than that of V-30" . 

The defect concentrations as a function of temper- 
ature and composition x were calculated with the in- 
dependent clusters approximation. [3, H3] This method 
is exact at low temperatures low defect concentrations 
where atomic vibrational effects and correlation among 
defect clusters are negligible. Under the constraint of 
x = 2([Vu] - [Iu]) + [Io] + Ei nPi + 2 Ej ",l>. ~ Wo] 
where i runs over COT-v and COT-o clusters |20l| and j 
over V-30 and V-40 (with n is the excess number of 
oxygen in each cluster, and the coefficient 2 before the 
second summation arises from the fact that these clusters 
are defined on the oxygen sublattice), we then got the 
defect concentrations by Eq.(l) of Ref. [20| in the closed 
regime where no particle-exchange with the exterior oc- 
curs. 

The results at 1500 K are shown in figure [2j where all 
unshown defect concentrations are small enough to be 
ignored. Cluster V-30 and V-40 compete each other 
throughout the whole x region, while their concentra- 
tions are always 1 order smaller than uranium vacancy 
and 2 orders smaller than that of Oi and COT-o. At 
lower temperatures their concentrations reduce further, 
and can be neglected. 

In this way a physical picture about oxygen defects 
appears, i.e., the nontrivial oxygen interstitials in hyper- 
stoichiometric regime are only Oi and COT-o. All O 
and O should come from COT-o clusters. This picture 
describes very well the high x region where U4O9/U3O7 
dominated. [20| It also explains the origin of asymmet- 
ric O and O interstitials and the non-negligible O v at 
T^ 300 K and x ~ 0.1 qualitatively. However, it left 
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a quantitative discrepancy: the experimental occupation 
number of O' and O" is 0.08:0. 16@ (which had been 
corrected to 0.13:0.12 late 3]), 0.14:0.12, and 0.33:0.10[| 
when x at around 0.11—0.13. It evidently asks for a low 
O :0 ratio, but COT-o cluster can provide only a ratio 
of 12:1. [13 

This disagreement is unlikely due to the failure of 
the current model. Three independent measurements on 
U4O9 and U3O7 has definitely proved that the clusters in 
these phases are COT exclusively, i.e., the O :0 ratio 
should approach to 12:1 at high x region. Meanwhile, 
the monotonicity of defect concentrations vs x is insensi- 
tive to the exact formation energies. It seems impossible 
that V-30 /V-40 cluster (contributor to O ) would 
have a considerable concentration at x ~ 0.11 but disap- 
pear at a higher one. However, the possible decomposi- 
tion/reassembly of COT clusters might have some impact 
on this issue, and we look forward to some understanding 
from this aspect. 

On the experiment side, ambiguity existed in its post 
data analysis. Recall the correction of the O :0 ratio 
from 0.08:0.16 to O.13:O.120[| and the similar vagueness 
in U4O9/U3O7 about the position of the center oxygen 
in COT clusters, [1, ©, Q3, ll| we know it is very diffi- 
cult to extract the explicit atomic geometry from powder 
neutron diffraction measurements, especially for VO2+X 
that is not an ordered phase and some peaks were broad- 
ened drastically. Another cause that might affect the 
O occupation is oxygen partial pressure. If experiment 
were performed in an open atmosphere and the powder 
were fine enough, then the above model need to be re- 
vised to take oxygen-exchange with the exterior into ac- 
count. Then the defect concentration becomes propor- 
tion to[IH pi ~ exp^-T^) (Po 2 f(T)) rH/2 where oxy- 
gen partial pressure Pq 2 could modify the concentration 
of V-iO /V-40 cluster greatly. Anyway, it requires 
more endeavor from theorist and experimentalist to re- 
move this discrepancy. 

In summary we clarified the defect clustering struc- 
ture in UO2+Z by first principles calculations. The Willis 
clusters have been proved generally failed. A picture of 



oxygen defective UO2 was established based on thermo- 
dynamical competition between O^ and COT-o cluster, 
which explains most available experiments very well. The 
last confusion is about the role of O^ in U02+ x - Willis 
once argued that oxygen should not occupy this site 
based on his observations, [l], [|J Above discussion indi- 
cated this argument is inappropriate. In order to explain 
the Willis' observation that Oi is negligible we should 
take the effect of temperature and composition into ac- 
count. Using the calculated defect concentrations, we 
plotted a pseudo phase diagram for oxygen defects in 
UO2+1C as shown in figure[31 Note each COT-o cluster 
contributes 13 single oxygen interstitials and the figure 
was plotted according to which kind of interstitial is ma- 
jority (from Oi or from COT-o). The figure clearly il- 
lustrated the predominant region of O^ is at low x and 
high temperatures. Mind all experiments available so 
far were conducted in the region marked by the hatched 
area, [2j, y, |4J which is far from the O^ territory. It is then 
understandable why experiments failed to detect Oj. 
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FIG. 3: Pseudo phase diagram for oxygen interstitials in 
UO2+1, where the hatched area marks the region in which 
the neutron diffraction measurements performed. 
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